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Abstract 

We present a theorem in d-dimensional static, spherically symmetric spacetime in generic Love- 
lock gravity coupled with a non-linear electrodynamic source to generate solutions. The theorem 
states that irrespective of the order of the Lovelock gravity and non-linear Maxwell (NLM) La- 
grangian, for the pure electric field case the NLM equations are satisfied by virtue of the Einstein- 
Lovelock equations. Applications of the theorem, specifically to the study of black hole solutions in 
Chern-Simons (CS) theory is given. Radiating version of the theorem has been considered, which 
generalizes the Bonnor-Vaidya (BV) metric to the Lovelock gravity with a NLM field as a radiating 
source. We consider also the radiating power - Maxwell source ( i.e. {F^yF'^'^Y , q = finely - tuned 
constant ) within the context of Lovelock gravity. 
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I. INTRODUCTION 



The string theory motivated higher dimensional gravity, known as the Lovelock gravity 
|l| attracted much interest in recent years This theory is known to admit the most 
general higher order invariants in such combinations that the field equations preserve their 
second order form. These are important features concerning divergences at smaller scales 
and ghost free structure toward a quantum theory of gravity. The first few orders of the 
Lovelock Lagrangian are well known: zeroth /first order is just the cosmological constant 
(A) /Einstein-Hilbert (EH) term. The second order Lovelock term is known also as the 
Gauss-Bonnet (GB) term, which consists of quadratic invariants and is a highly non-trivial 
contribution, especially in higher dimensions js]. Going to higher order corrections serves only 
to add contributions at higher levels of intricacy which are crucial in a possible quantum 
gravity. In spite of all complications it is remarkable that in a spherically symmetric line 
element exact black hole solutions are found in higher order Lovelock gravity [3]. These 
metrics are sourced by Maxwell and Yang-Mills (YM) fields. More recently, we obtained 
black holes in a power- YM source in Lovelock theory j^. By the power- YM source it is meant 
that the source consists of a power of the YM invariant, i.e., [F^^F°'^'^Y . Here F^^ refers 
to the YM field with gauge group index a and g > is an arbitrary constant parameter. 
Although q may stand arbitrary, energy and causality conditions restrict it to a certain set 
of admissible parameters Sj. 

By the same token, in this paper we consider Lovelock gravity of all higher orders coupled 
with a non-linear Maxwell (NLM) source j^. We prove a Theorem, covering Lovelock terms 
to all orders - albeit with proper constant coefficients - which generates new metrics once 
the NLM energy-momentum is known. The energy-momentum is given by C (J-") , in which 
J-" = F^^F'^" denotes the Maxwell invariant. The field equation satisfied by F^,^ will be 
referred to as the NLM equation. The Theorem involves spherically symmetric metric ansatz 
together with the NLM Lagrangian which leads to a general class of metrics. Static electric 
fields fall within the range of our Theorem which yields Born-Infeld (BI) electrodynamics 
as a particular example. A static, pure magnetic YM field constitutes another example of 
BI type. 

For the first/second order Lovelock terms the Theorem addresses to EH/GB gravities, 
resulting in metrics expressed in solutions of first/second order algebraic equations. With 
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equal ease, we extend this result to any higher order Lovelock gravity in terms of a higher 
order algebraic equations As an interesting example we obtain Chern-Simons-Born-Infeld 
(CSBI) black hole solutions in odd dimensions. We elaborate on d = 5, to show in particular 
that the thermodynamically well-behaving CS black hole preserves, its good features when 
coupled with a BI source in the general relativity limit. We note that combination of 
CS and BI types is in an artificial manner since CS/BI black holes arise naturally in odd 
/even dimensions. Such a combination works only when the *-F^i/ (i.e., dual of F^^) vanishes, 
which occurs in a restricted type of sources. BI black holes are known to arise naturally from 
Pfaffians only in even dimensional spacetimes j^. Depending on the topological parameter 
X = 0, ±1 and A ^ we investigate the Hawking temperature of the 5-dimensional CSBI 
black holes. Finally, we investigate the implications of our Theorem in the Eddington-Bondi 
form of radiating metrics Q]. This is the time dependent version of both the mass and charge 
so that the metrics become time dependent. That is, we extend the Bonnor-Vaidya 8|| form 
of radiating metric to the d-dimensional Lovelock gravity coupled with a NLM field. We 
consider various types of non-linearities in electromagnetic field as particular examples. 

Organization of the paper is as follows. The Theorem, its proof and applications are given 
in Sec. II. Generalization of the Theorem to the radiating metrics in higher dimensions with 
NLM as source is studied in Sec. III. In section IV, we investigate the case of power-Maxwell 
non-linearity. The paper is completed with the conclusion in Sec. V. 



II. A THEOREM FOR SOLVING EINSTEIN-LOVELOCK-NLM EQUATIONS 

Theorem : Let the d-dimensional static spherically symmetric spacetime sourced by 
non-linear electromagnetic field be described by the action {8nG = 1) 

S=^j 6/x^v^|- ^'^~^^^^'^~^^ A + £i + a2>C2 + ... + a[^]i:[<^] +/:(^)| , (1) 
in which 

£„ = 2-"C'r"rf'"^^''^ft ■■■R'"n\ , ^ > 1, (2) 

" C\...Cna\...an CL\0\ CLnOn^ — ' V / 

is the nth order Lovelock Lagrangian, an is a real constant and the bracket [.] refers to 
integer part. Herein T = F^^F^^'^ , is the Maxwell invariant for a static field 2-form 

F = E{r)dtA dr. (3) 
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If the C{JF) satisfies the non-linear Maxwell (NLM) equation 



with Cjr = ^ and *F the dual of F, then the Einstein equations admit the solution 



ds^^-{x-r''H {r))dt^ + 



-dr^ + r^dnl_2 



(4) 



(5) 



where x = 0, ±1 and in which H (r) is the solution (or solutions) of the following algebraic 
equation of order [^^] 



E {rf = - + 



fe=i 



3 r'^-i {d - 2) r"^- 



jjr"-' {jr-2jr^T)dr. 



(6) 



2k 



Here m is an integration constant, au = 11 ("^ ~ and ai = 1. 

j=3 

Proof: Variation of the action with respect to the metric tensor g^i, yields the field 
equations in the form 



\ 



\ 



k=2 

where Qj^^^^^ is the Einstein tensor, while 



(d-2)(<i-l) 



(7) 



A* 2.^+^ bci...Ckdi...dk aibi'" a^bk' 



_^aai...akbi...bk j^cidi j^Cf^df, 



The T^'' is given by 



(8) 



(9) 



which clearly gives T^* = T^*" = |£ — C-jT ^ stating that = Q/ and T^^' = Now we 
introduce our metric as given by (5), where the choice of gu = — {grr)~^ is a direct result of 
= Q/ up to a constant coefficient where we set it to be one. By starting with the line 
element (5) one gets 



E«fc {r'-'H{r) 

k=2 



and 



^ ('■"'f^ w)" - ^ E 3' (--"^ w')" . 

k=2 



(10) 



(11) 
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where (.)' = f (■) ■ Eq.s (10) and (11) admit 

-r'-'g'% (12) 



dr 

„d-2 



d-2 

and imposing Einstein equations yield 



d-2 \ ' 6 
This is equivalent to, 



/I \\' /I 



-C-C^J- =r'-' , (14) 



which integrates to 



d-2\2 \2 



r'^~^LT^\T\ = constant. (15) 



This result is in conform with the solution of the NLM equation (4) and the integration 
constant can be identified with the electric charge. Now recall from (7) , that 



, , {d-2){d-l) 



Gt+- ^A = 7;*, (16) 



or equivalently 



fe=2 

(17) 

which implies 

KM 1 

5: a,E {rf = 3 + - (^3^/^'" ' ^^^^) (18) 

for the integration constant M = j^^)- ^^^^ completes the proof of our Theorem. 

Remcirk 1 In the Theorem, we considered only a static electric field whose 2-form is F — 
E (r) dt A dr, giving T — — 2E'^ . Since T is only a function of r so is the Lagrangian C [F) , 
this aided in the proof of the Theorem. We must add that the Theorem becomes of practical 
use whenever the expression on the right hand side is integrable. Otherwise we should resort 
to the multipole expansions for the asymptotically flat spacetimes. 
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Remark 2 In the case of pure magnetic field the Theorem is applicable only if J-" = F^^F^^ 
is only a function of r. This implies 

d-2' 

which leads to the same metric function in the form 



tA =r'~V% (19) 



AM 2 

k=l 

Remark 3 In the case of a general energy - momentum tensor 



Y: a,H (r)' = ^ + ^ - 'dr. (20) 



T;^diag(T,',T/,T,;\...), (21) 

in which 

T/ = T/,r,f^ = r,f^ = ... (22) 



and 

^— ^T/j =r'^-%f% (23) 
again, a solution in the form of (5) satisfies the Einstein equations with H (r) given by (20). 

Remark 4 The case of C — ICjT = must be excluded, since it implies a Lagrangian of 



the form C =\/J-', which fails to satisfy the energy and causality conditions fJ/. This form 
of the Lagrangian lacks also the linear Maxwell limit. 

Example 1: As an application we consider the case of pure electric Einstein-Born- Infeld 
(EBI) black hole solution. The pure electric BI Lagrangian can be written as 



/:(^)=4/3^|l- Jl + ^j, (24) 

where J-" is the electric field invariant given by = 2Ft^F*^ = -2E {ry . Note that for 
/3 — 7- oo we recover the standard, linear Maxwell Lagrangian. Consequently 



Cr = - (25) 



1 + ^ 



and therefore upon solving the non-linear Maxwell equation one finds 



Finally we obtain 



3^7^- JI^^^5r'-'% 'ir == - + (^3^- (27) 



(d - 1) (d - 2) I y /32^2(d-2) y {^d-l){d- 3) r2(rf-2) 
^/l d-3 3d -7 g2 



,2' 2(c/-2)' 2(rf-2)' /32^2(d-2)y' 
in which stands for the hypergeometric function. 

Example 2: Another example for the case of non-electric field is given by the Einstein- 



Yang-Mills (EYM) non-linear electrodynamics black hole solution. In fact in Ref. jo] we 
find that 

V.T,^2,-'(.^f:^I^^l (28, 

and 



j ^4 A ^ (d-2K^^^3)Q^ ' ^ ^ 

which clearly satisfies the conditions of the Theorem and the Einstein equations admit a 

black hole solution with 

, , A 4m 4/32 

H (r) = 1 h 

^ ' 3 (d-2)r'^-i (rf-l)(rf-2) 



(.-2).^-T^-^-^"V^^+ 2/32 ' 

in conform with the solution given in Ref. {q]. 

Example 3: Our next example will be the general form of energy momentum tensor given 



by Salgado 



10| which states that 



C 

KiDian.) = „n-k) [1)1)^)-)^]) {C,k: constants), 



-fj.{Diag.) ^n(l-fc) 

admits solutions for Einstein equations. Now we show that this is a natural result for the 
case of Rem. 3. In Rem. 3 let's consider 

T/=T/, T,f = A;T/, (31) 

then, from equation (23) we obtain 

'-tA =kr^-%\ (32) 



d-2' 



or, in a straightforward calculation one finds 



C 



j,n{l—k) '■ 



(33) 



where C is an integration constant. This verifies that the Theorem proved by Salgado 10 | 
turns out to be a particular case of our more general Theorem. 

Example 4 '■ One may notice that a proper choice of ak leads to a Chern-Simons (CS) 



ll| gravity in odd dimensions. To do so we set 



ttfe = — , for /c > / and — — = — , 

«! 3 «! 



(34) 



and we rewrite (18) as 



fc=0 

Now we choose 



where 



^ 2 J /Ml \ 



A; 



3 a, [^]' 



to get from the binomial expansion 



+ (r))l 
This implies that 

H (r) = + — 

y ^ £2 ^ £2 



M 



A ao 



M 



r'^-i (d - 2) r' 



(c/ - 2) r'^- 



(35) 

(36) 
(37) 

(38) 
(39) 



where a = +1 if [^2^] is an odd integer and o" = ±1 if [^2^] is an even integer. 
The latter equation for d =odd, x = 1 admits (£^ > 0) 



1 M 



2B (r) 



2 

d-l 



(40) 



(rf-2), 

where / (r) = (x — r^H (r)) is the metric function and B (r) = f^z'^~'^Tf (z) dz. This metric 
function by using (39) and (40) becomes 



m + 1 



{d-l)B{r) 



2 

d-l 



(41) 



(c/-2)£('^-3)_ 

in which the new integration constant m is related to M [H]. In the sequel we investigate 
some thermodynamic properties of this solution. 
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1. With a = +1 or is an odd integer (d = 7, 11, 15, ...) 



In this case one finds the Hawking's temperature as 



H 



+ 



{d-2) (£2 + r2)l — 
The specific heat capacity of the black hole for constant charge is defined by 

dS 



H 



and is obtained as 



where 



4(¥)! 

(d-l) (ci-2)7r('^)rt'T 



T = (d - 2) + 



* = (d-2) + 



d 



r+ {f + rl) -^Tl (r+) + T* (r+) (3f - (d - 6)) 



(42) 



(43) 
(44) 

(45) 
(46) 



2. With a = ±1 or is an even integer (d = 5, 9, 13, ...) 

It is clear that in this case a — —1 does not claim any horizon and therefore it is out of 
our interest but for cr = 1 branch the Hawking temperature and the specific heat capacity 
are given as (42) and (44) but there exists an additional constraint on the free parameter in 
order to have a real metric function. 

To complete this section we give an example for the CSBl black hole (x = 1, cr = 1) in 
5-dimensions. To do so we recall that in 5-dimensions 



B (r) = fz%' {z) dz 



1 + 



3g2 



^2^6 y + 4 ^2 2^1 1 ^ ' o ' o ' 



1 1 4 



2' 3' 3' 



(47) 



and therefore 



2^4 



2/3'r 



1-1 1 + 



^2^6 



r ^ /'I 1 4 



2' 3' 3' 



(48) 

In Fig. 1 we plot / (r) for the specific choices of mass, charge and the cosmological constant. 
The corresponding temperature and specific heat capacity plots are given in Fig.s 2 and 3. 
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Example 5 (Clouds of strings as source): We consider another application of the Theorem 
;hat incorporates energy- momentum tensor representing clouds of string type matter fields 
I2I. More recently [l3|, clouds of string type energy- momentum is considered in Einstein- 
Gauss-Bonnet (EGB) gravity. In Rem. 2 and 3, we considered the case for a general energy 
- momentum tensor and its corresponding solution. Clouds of string type matter fields obey 
the condition imposed on the energy - momentum tensor stated in Remark 3, such that, 



(49) 



real constant, 



which leads 



Tl: = 0. (50) 
Therefore the general solution for the metric function is obtained from the algebraic equation 

2a 



2^ (^kH (r) = - + 



k=l 



(51) 



The solution for H (r) generalizes the solution obtained in Ref. [3| to higher order Lovelock 
gravity. And hence, our general solution includes the solution obtained in 13| if we restrict 
the spacetime dimension to d = 5. For this particular case the solution is obtained from 



H{r) + a,H{r)' = - + --— , 



or equivalently 



/ (r) = X + 



1± Wl 



A M 



2a 



"2 - + -r - TTT 



(52) 



(53) 



This solution is nothing 
cloud type matter fields 



4^2 \ V V 3 r4 3r^ 

)ut the black hole in EGB gravity in the presence of the string 



12| . In this theory one can easily check that the corresponding 



energy momentum tensor in d-dimensions is given by 

a a 



diag 



(54) 



At this stage we wish to go a step further to relate our solution to the CS solution in odd 
d- dimensions as (x = 1), 



fir) 



r 



a 



m + 1 — 



{d-2) 



-3) 



2 

d-l 



(55) 
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which in 5- dimensions becomes 



/W = 1 + ^±V"^ + 1-^- (56) 

The thermodynamic properties of this solution is also investigated. The Hawking tem- 
perature ( Tu) and heat capacity Ca at constant string parameter a is calculated at the 
location of the event horizon ( ) are given by, 

_ /2(3^^ + a) + 3r^ 

- 37r/2(/2 + rD ' (^^) 

^ 3^rl{P + rl)r {Sr, + a)+3rl] 

2[P{Qrl + 3P -2arh) + 3rl] ^ ' 

We also analyzed the thermodynamic stability which is indicated by the positive heat 
capacity Ca- If the heat capacity has unbounded discontinuity at particular points of r/j, 
this implies possible phase change from stable to unstable black hole solution. As illustrated 
in Fig. 4 the transitions from stable to unstable black hole solution is not continuous and 



therefore possible Hawking-Page type phase transition occurs [15|. The occurrence of the 
phase transition crucially depends on the ratio of a/ 1. In Fig. 4, the values of this ratio that 
creates phase transitions is depicted. 

III. A GENERALIZATION TO THE RADIATING METRICS 

Following Bonnor and Vaidya (BV) js], we consider a non-linear electrodynamic La- 
grangian C (J-"), a null current J and a coupling term A^J^ added to the original Lagrangian 
such that 

(i(*F£^) = *J. (59) 
Now, we consider the d-dimensional version of the BV metric 

ds"^ = -{x- r^H (r, u))du^ + 2edrdu + r^dQl_^, (60) 

with the outgoing null coordinate u and field 2-form 

F =E{u,r)dr Adu. (61) 

This gives 

*F =E{u,r)y/^d9iAd92... Ad9d-2, (62) 
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in which g — det {gnu) and the Einstein equation is given by 



(63) 



Here 



and 



k=2 

rpv{fluid) _ _YVY 



(64) 



(65) 



for a null vector V^. We start now with NLM equation (50) which leads to 

d (*F£^) = d{E {u, r) C^^^dOi A d^2-- A ddd-2) = 
[{E {u, r) £TV^)r dr + {E (u, r) Lr^f^)^ du\ A dQ^ A dQ^- A dQe.-i J- 

By using the relation between 1-form current J and its dual i.e., J = (— 1)'^ ** J one finds 



(66) 



J= (-1)^ {{E {u, r) Lt^), * {dr A dO^ A d^s-. A d9d-2) + 
{E {u, r) Ljr^f^)^ {du A ddx A dd^... A ddd.-2)\ ■ 



Prom the metric we find 

^{dr ^d^x^d^2...^dQA-2) 



■1) 



'-9 



{dr — efdu) , 



(-1) 



d-i 



-du, 



and therefore 



{du A d9i A d92... A d6'd_2) = 



J= {E {u, r) C^r^-')^ ^ {dr - efdu) - {E {u, r) £^)„ du. 



(67) 

(68) 
(69) 

(70) 



This current is going to be null i.e., — (7^,0, ...,0) = J„(5^ which means that 
{E (m, r) Cjrr^~'^)^ = integrates to 



E{u,r)jC^r''-^^Q{u) 



or 



(71) 



(72) 



12 



for a u dependent charge Q{u). After considering these results we find 

Q{u) 



F = 



-dr A du, 



3=-(-l 



T 



{E{u,T)Lj:)du = {-\) 



d-2 



du. 



where Q (u) = The exphcit form of the Maxwell field can be expressed by 



= 2Fru (g'^'g^^'F^p) 



while our metric tensor g^j^^, and g'^'^ are 
/-/ 6 . 



g|^u 



\ 



e 

r^sin^^ . 



V 



I 






1 

e 








1 

e 


/ 














1 














1 

sin^ 6 



\ 



giving 



The energy - momentum tensor components are given by 

rpr(em) i 



rpOi{em) 1 ^ 

rpr{em) rpu{em) g 



The null- fluid current vector is V^= {Vu, 0, 0) = VuS^ and therefore 
implying that 

g^.V^V = grr {V'f , 

which obviously vanishes and 

rp.{fiuid) ^ ^ _vrvj';5^ = -e (K)' S-^S^. 



u{em) 
u , 
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Finally, we give the explicit form of the energy momentum tensor as 

( \{C-2CrJ^) . A 

-e{Vuf \{L-2LtT) 0.. 



rpv 














. . 
i£ . . 



V 

Similarly the Einstein tensor can be expressed by 



(85) 



d-2 



\ 



1=1 

m 







d-2 I ^d-l 



i=l 













E 



i=l 



where 



/ (m, r) = 1 — r^i? (m, r) . 



(86) 



(87) 



After all these arrangements we are ready now to solve the Einstein equations. We start 
with the uu component 



1=1 



i=l 



{d - 2) r" 



(88) 



One can easily check that the rr component gives the same result while the angular compo- 
nents give 



r'^-^ E c^iH' I = ->C ^ ( 3-— {C - 2Cj.J^) ) = r'^-^£ (89) 



i=l 



d-2 
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which is nothing but the NLM equation (14), already satisfied. In the last step we work on 
= T!" giving 

(] -2 8 

— ^5^ES,ff = -.(l4f (90) 

and therefore 



d-2jM{u) 2 r J-2^r 



2 ■ „_2),.-. (91) 

which finally determines the null-fiuid current vector by 



IV. RADIATING POWER-MAXWELL SOURCE IN LOVELOCK GRAVITY 

In this section our choice for the NLM field consists of a particular kind, namely 

C (J^) = -J^i (93) 
in which q stands for a constant parameter and J-" =F^ijF^^ , as before. Let us note that 



3- 



this particular version of non-linearity attracted considerable interest in recent years 
Although the parameter q (i.e. the power) may be assumed arbitrary, imposition of the 
energy - conditions with other requirements restrict g to a limited set of values. The rest of 
the gravitational action will be chosen as in the previous sections. Similarly, our choice of 
the line element follows that of (60). The electromagnetic field 2- form turns out to be 

F = -^^dr A du, (94) 
and energy-momentum component ^ 

T/ = ^(/:-2£^^)= (95) 



so that 



g2^2(d-2) 

Following the procedure of the previous section further, which led us to the condition 
(88) we obtain now, in analogy 
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[ d-l 
2 



J: a,H' = ^ - hAu)/-^ (97) 

i=l f 



in which 



W»)=(-l)'(^)(^) ■ (98) 



'2g - 1\ f2Q'^{u)\ 29-1 

It is observed that an arbitrary q does not guarantee the reahty of the metric function, which 
enforces us to choose q appropriately. 

The particular dimensionality d — h brings in significant simplicity to the foregoing 
expressions such as 

aiH + a2H^ = - h5{u)r'-^ , (99) 

where ai — 1 and 0:2 = 2a. The quadratic equation for H can easily be solved and the 
corresponding metric function /(r, u) is determined by 



f{r, u) = x + ^U±Jl + Sa - h{uy-^ ] ] . (100) 



In analogy to Eq. (82), the null-current component for the present case turns out to be 



in which, as usual, a 'dot' represents 

To be able to proceed further with the thermodynamical properties in the present choice 
of d = 5 ( and x — +!)> we must determine the apparent horizon through /(r, u) — 0. This 
leads us to the algebraic relation 

2a- M{u) + rl + h5{u)r'l = (102) 
in which rh denotes apparent horizon (if any) and we have abbreviated 

k^2('^). (103) 



.2g- 1, 

Tab. 1 shows the relation between q and k for certain leading numbers of interest. With 
reference to Table 1, we can find a finely-tuned set of q powers so that h5{u) will be real. 
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Clearly, any odd/even integer q will do the job whereas non-integer g's will not serve the 
purpose. For instance, g = 4 is a good choice. Accordingly, we find 

^(0!M)'%o, ,104, 

and 



{u) = ^ (Vl + Ah^{u) {M{u) - 2a) - l) (105) 
M{u) > 2a. 



The corresponding metric function takes the form 



/(r, n) = 1 + ^ ( 1 ± a/ 1 + 8a^ - h,{u) 1 . (106) 



I 1± A/l + 8a— , 
4« \^ V r4 

Different choices of q values from Table 1 can be treated in a similar manner to obtain 
the corresponding f{r,u) function, which we shall not go any further in this paper. 



V. CONCLUSION 



Methods to generate new solutions for various energy-momenta are available in the liter- 
ature. In this regard, recently we have proved a theorem that generalizes Salgado's theorem 

. The whole issue in such a problem is the physical significance of the energy-momentum 
under consideration. Scalar field source, for instance, is known to work only in d = 4. Lin- 
ear Maxwell electrodynamics has already been incorporated in higher dimensional Lovelock 
gravity . Our example 5 which is a particular version of the Theorem concerns a cloud 
of strings without Maxwell field. We explore thermodynamic properties of such a cloud for 
a particular solution in d = 5 spacetime dimension, representing Chern-Simon black holes 
which undergoes a Hawking-Page phase transition (Fig. 4). In this paper we extended the 
problem to cover a non-linear Maxwell (NLM) source. Among this class of theories Born- 
Infeld (BI) is the most familiar one, which arises as a particular example to the Theorem. 
The Chern-Simons-Born-Infeld (CSBI) metrics that we found by using our Theorem turn out 
to be thermodynamically stable. Next, we address to dynamical problem that covers time 
dependent metrics and NLM fields. Stated otherwise, we have generalized the well-known 
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dynamical Bonnor-Vaidya (BV) metric to the general Lovelock theory with NLM sources. 
A radiating null current source naturally accompanies the radiating energy - momentum of 
such metrics which lose mass and charge. It is needless to remark, finally, that the opposite 
problem of the 'shining star', namely, the collapse of time dependent energy - momentum 
and null radiation current is also solved in the same theory. 
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Captions: 

Table 1: The relation between the parameters k and q according to Eq. (103). 

Figure 1: The plot of / (r) (Eq. (48)), for the specific parameters m = £ = q = l, as/3 
ranges from to oo. It is seen that black holes solutions are available only for (3 < ^critical = 
0.2276. 

Figure 2: Hawking temperature Th versus event horizon radius r+ for the specific param- 
eters £ — q — 1 and for /3 e (0, ^critical] ■ By taking the absolute value of Th automatically 
deletes the negative temperatures as non-physical. 
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Figure 3: specific heat capacity Cq versus r+ for versus r+ for £ = g = 1 and ^ e 
^critical] ■ This plot, (together with Fig. 1), reveals that our CSBI black hole solution is 
thermodynamically stable. 

Figure 4: Heat Capacity Ca versus the horizon radius for specific parameters a = 5 
and £ = 1. The singularity in Ca and therefore occurrence of Hawking-Page phase transition 
is clearly seen. The dash region in the inscribed figure depicts | versus for which such a 
transition occurs. 
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